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Abstract 

Two counter-propagating cool and equally dense electron beams are modelled with particle-in- 
cell (PIC) simulations. The electron beam filamentation instability is examined in one spatial 
dimension, which is an approximation for a quasi-planar filament boundary. It is confirmed, that 
the force on the electrons imposed by the electrostatic field, which develops during the nonlinear 
stage of the instability, oscillates around a mean value that equals the magnetic pressure gradient 
force. The forces acting on the electrons due to the electrostatic and the magnetic field have a 
similar strength. The electrostatic field reduces the confining force close to the stable equilibrium of 
each filament and increases it farther away, limiting the peak density. The confining time-averaged 
total potential permits an overlap of current filaments with an opposite flow direction. 

PACS numbers: 52.38.Hb,52.35.Qz,52.65.Rr 
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The electron beam filamentation instability (FI) generates magnetic fields in energetic 
astrophysical^' ^ and solar flare plasmas^ and in laser plasma interactions^'-, if the beam 
speeds |vb| are comparable to c and if the densities of the counterstreaming beams are 
similar^. It has been investigated with one-dimensional (ID) PIC and Vlasov simulations^' - 
and with two-dimensional (2D) PIC simulations*^'-—. The counterstreaming electron beam 
instability has also been examined with a 3D PIC simulation^^ . Mobile ions and a guiding 
magnetic fleld have been taken into account^'—'— and statistical properties of the FI have 
been obtainec^^'— '— . 

The FI triggers the growth of waves with the wavevectors k ± v;, over a wide band 
of A; = |k|, where the wavenumbers k are of the order of the inverse electron skin depth. 
The electrons are deflected by the magnetic fleld perturbation, and electrons moving in 
opposite directions separate in space. The net current of these flow channels amplifles the 
initial perturbation and, thus, the tendency to form current channels. The magnetic fleld 
amplitude grows exponentially and it saturates by the magnetic trapping of electrons^. The 
FI can also couple nonlinearly to electrostatic waves^'— '— . It has been suggested^'— that 
it is the magnetic pressure gradient that gives rise to the electrostatic fleld that grows, when 
the FI saturates, but it has not yet been demonstrated quantitatively. This is the purpose 
of this paper. 

We consider here the FI driven by equally dense and warm electron beams, which have 
a Maxwellian velocity distribution in their rest frame. This case is important, because the 
growth rate of the FI is highest relative to the competing mixed mode and two-stream 
instabilities for symmetric beamsS. We study the FI with a particle-in-cell simulation code^ 
that is based on the electromagnetic and relativistic virtual particle-mesh methoc^i^. 

The FI is modelled in a simulation reference frame, in which both beams move into 
opposite directions at the speed modulus Vb = 0.3c. We isolate the FI by selecting a ID 
simulation box that is oriented orthogonally to the beam velocity vector v;, and we resolve 
all velocity components. This is an approximation for a quasi-planar boundary between 
fllaments with an oppositely directed electron flow. They occur in warm plasmas, if the 
conflning magnetic fleld cannot overcome the thermal pressure and they are characterized 
by planar magnetic flelds^'— '— . The periodic boundary conditions of the short simulation 
box result in the development of only one pair of fllaments. The restriction to one dimension 
inhibits the merging of the fllaments^ and we can analyse the relation between the electric 
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and magnetic fields of tlie quasi- stationary filaments. 

Beam 1 has the mean speed v^i = VbZ and the beam 2 has Vh2 = —VbZ. Both beams 
are spatially uniform and have a Maxwellian velocity distribution in their respective rest 
frame with a thermal speed Vth = {kbT/rrif.)^'^ of Vb/vth = 18. The ID simulation box 
with its periodic boundary conditions is aligned with the a;-direction. We thus denote 
positions by the scalar x. The plasma frequency of each beam with the number density 
rie is Up = {e^rie/rrieeof'^ . The total plasma frequency Vtp = ^/2ujp. The electric and 
magnetic fields are normalized to Ejsi = eE / cmStp and = eB/mS^p and the current 
to Jn = J jlueec. The physical position and time are normalized as xat = x/A^ with the 
electron skin depth = cjVtp and ^at = tfip. We drop the indices N and x, t, E, B, J are 
specified in normalized units. 

The box length L = 0.89 is resolved by Ng = 500 grid cells with the length A^.. The 
simulation time ts = 125. The phase space distributions fi{x,p) of beam 1 and f2{x,p) of 
beam 2 are each sampled by Np = 6.05 ■ 10'' computational particles (CPs). The total phase 
space density is defined as f{x,p) = fi{x,p) + f2{x,p). 

The electrons and their micro-currents are redistributed by the FI along x. The charge- 
and current-neutral plasma is transformed into one with Jz{x) ^ 0. The z-component of 
Ampere's law is in the ID geometry d^By = Jz + dtE^. A oc sin {kx) gives a (x sin (kx) 
and By oc — cos (kx) so that and By will have a phase shift of 90°. Figure [1] reveals this 
phase shift between By and the evanescent Ez- It also shows, that an electrostatic -Ea;-field 
grows. The By{x,t) and the Ex{x,t) oscillate in space with the wavenumbers ki and k2, 
respectively, where kj = 27Tj/L. Both fields are spatially correlated. The comparison of Ex 
and By aX t = 56 demonstrates that E'^. = if -Bj^ = or if d^By = 0. 

We determine now the relation between E^ and By. Let EB{x,t) be an electric field 
along X, which excerts the same force on an electron as the magnetic pressure gradient 
force does. This electric field is given in our normalization (charge q = —1) as EB{x,t) = 
—By{x,t)dxBy{x,t). We note that By{x,t > 56) is quasi-stationary, while Ex{x,t > 56) 
oscillates in time. The oscillation amplitude of Ex{x,t > 56) is approximately constant 
and it apparently oscillates around a stationary background field. It is helpful to average 
the Ex{x,t) and the EB{x,t) over the time interval ti = 56 to t2 = 125 to give Ex{x) = 
{t2 -tiY^ j'll Exix, t) dt and Eeix) = (ts - ti)"' j'' Eb{x, t) dt. 

Figure [2]( a) displays the Ex{x,t = 56) when the FI has just saturated and reached its 
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FIG. 1: (Color online) (a)-(c) show By, and E^. All fields are stationary in space and By is 
quasi-stationary also in time for t > 56. E^ oscillates in space twice as fast as By and both are 
spatially correlated for 55 < t < 125. The phase of E^ is shifted by 90° relative to that of By and 
it is damped. The By, Ex fields at t = 56 are compared in (d). 
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FIG. 2: (a) The Ex{x,t = 56) (dashed curve) and 2EB{x,t = 56) (solid curve), (b) The time- 
averaged Ex (dashed curve) and Eb (solid curve), (c) The number densities for t = 56, normalized 
to 2ne, of both beams separately (beam 1 is almost confined to 0.2 < x < 0.6) and both densities 
added together. 

peak amplitude and it compares it with EB{x,t = 56). It turns out that Ex{x,t = 56) ~ 
2EB{x,t = 56). The time-averaged fields fulfill Ex{x) ^ Eb{x) in Fig. EJ^b). The Ex{x,t > 
56) oscillates in time with an amplitude ~ Eb{x) around a stationary background field 
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FIG. 3: (Color online) The 10-logaritlimic phase space densities in units of CPs at t = 56: Panel 
(a) shows the f{x,pz) with pQ = meVbT{vi)). The temperature and the mean velocity along z of 
the electrons are unchanged as a function of x. The density oscillates by the factor ~ 10. The 
fi{x,px) is shown in (b) and the f2ix,Px) in (c). 



with the amplitude ~ Eb{x). Both amplitudes add up to 2Eb{x) at t = 56. When the 
oscillatory and the background electric field have a phase shift of 180° in time, they result 
in a Ex{x, tc) ~ 0, for example when tc = 75 in Fig. [T](b). 

The Eb{x, t = 56) and Eb{x) correlate well in Fig. [2](c) with the normalized number den- 
sity distributions ^1^2(3^) = (2?^e)^^/ fi,2{x,p) dp of each beam and also with the summed 
distribution ni{x) + n2{x) at t = 56. The total density is modulated by about 30%, while 
that of rii and n2 varies by an order of magnitude. 

Figure [3] shows the electron phase space distributions f{x,Pz), fi{x,Px) and f2{x,Px) at 
the time t = 56. The mean velocity along z of the electrons of the beams 1 (i=l) and 2 
(i=2) is practically constant as a function of x. Any spatial modulation would be caused by 
the E X B-foTce, which is given by the product of E^ and By in our geometry. The effects 
of this force are small. 

The supplementary movie^ animates in time the evolution of fi{x,px) and fi{x,pz), 
where the color scale denotes the 10-logarithmic number of CPs. The electrons are redis- 
tributed along X but they keep their p^ almost unchanged. Their flow along x oscillates, 
giving a Jx{x, t) 7^ 0. The fi{x,Px) has a dense electron core, which rotates in the x,p3;-plane 
around x = 0.4. Two phase space vortices are convected with this rotating flow. The phase 
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FIG. 4: The fields Ej and the potentials Uj averaged over 68 < t < 125: (a) shows Ex (dashed), 
Ed = —vifBy (dash-dotted) and Ej- = E^+Ed (solid line). Positive Ej accelerate electrons into the 
negative x-direction. (b) shows the potential tlx (dashed), Ud (dash-dotted) and the Ut (solid). 
The potential at x = 0.4 is the reference potential. 

space motion of the electrons around the equilibrium points Xe with Ex{xe) = By{xe) = 0, for 
example Xg = 0.4 in the supplementary movie, reveals, that they are trapped by a potential. 

We can estimate the contributions of Ex and By to this potential after the satura- 
tion of the FI. We average the fields Ex{x) = {t2 — ti) ^ J^^ Ex{x,t) dt and By(x) = 
it2 -hY^ It' By{x, t) dt in t ime from ti = 68 to ^2 = 125. In what follows we consider 
beam 1 with Vb > 0. According to the supplementary movie, most of the electrons have 
the velocity components Vz ~ Vb and Vx <^ Vb- The electrons retain their initial Vy -C Vb, 
since no force component along y develops. The dominant component of the time-averaged 
magnetic force is thus Fx = VbBy for q = —1. The time-averaged force along x is then 
Fx = —{Ex + Ed) with Ed = —VbBy and we define Ft = Ex + Ed- The time-averaged 
potentials Uj{x) = Uqj + Ej{x)dx with the indices j = x,D,T are calculated from these 
fields and Uqj is set such that Uj{xe = 0.4) = 0. The potentials are given in Volts. 

Figure m displays the time- averaged fields and potentials. The Ex destabilizes the equiUb- 
rium position Xg = 0.4, because the negative Ex{x > Xe) close to Xe accelerates the electron 
in the positive direction and the positive Ex{x < Xe) close to Xg in the negative direction. 
The Ed{x ~ Xe) is confining the electrons around x ~ Xe- The \Ed\ > \Ex\ for x ~ Xg and 
Ft is thus a confining force. However, the electron acceleration at x ~ Xg is decreased by 
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FIG. 5: The trajectories of two selected CPs: (a),(c) show the x,px and x,pz diagrams of the CP 
1. (b),(d) show the corresponding diagrams for the CP 2. The circle denotes the starting point of 
the trajectory. Both CPs follow straight paths in the x^px plane for 0.33 < x < 0.47 and they are 
rapidly reflected outside this interval. 

Ex and increased at larger |x — Xe|. 

The CPs of the beam 1 should follow almost straight paths close to Xe = 0.4 and they 
should be rapidly reflected for |x — Xe| > 0.2. The potential difference Au = max{UT) — 

~ 1 /2 

min{UT) ~ 1700 V should trap electrons with speeds up to A^, = (2eAc//me) /fj, ~ 0.27. 
This matches the momentum spread of the cool core population in Fig. [3] and in the 
supplementary movie. The oscillations of Ex in Fig. [1] and, thus, of the strength of the 
confining potential explain the periodic release of electrons from this cool core seen in the 
supplementary movie. The oscillatory force imposed on the electrons by Ex{x, t) contributes 
to their heating. 

Figure [5] follows the trajectories of two CPs of beam 1. The circles denote the times, when 
the CPs start interacting with the fields and the trajectories are followed until t = 125. The 
CP 1 has a low initial modulus of Px and CP 2 a high one. Both CPs follow straight paths 
in the interval 0.33 < x < 0.47, in which Ej- in Fig. IHis small. The phase space path of the 
faster CP 2 is smoother than that of CP 1. The low speed of CP 1 implies a long crossing 
time of the interval with a low modulus of Et and the CP 1 experiences several oscillation 
cycles of Ex- Both CPs are reflected outside the interval 0.33 < x < 0.47 and they remain 
trapped, because they can not overcome the potential difference A[/. Both electrons change 
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their Vz only by a few percent and \vx\ < Vb/5, which is supporting our previous assumption 
of a dominant and constant force VbBy along x. 

In summary, we have examined the saturation of the filamentation instability (FI) driven 
by two counter-propagating, weakly relativistic and symmetric beams of electrons. The 
ID simulation box has been oriented orthogonally to Vf,. It can approximate the quasi- 
planar boundary between two filaments with oppositely directed flow, which shows up if the 
magnetic confinement cannot overcome the thermal pressure. This geometry is beneficial, 
because two of the three spatial derivatives in the Maxwell equations vanish, which separates 
the electrostatic and the electromagnetic fields. 

We have confirmed here, that the electrostatic field, which grows during the nonlinear 
phase of the FI and in a ID geometry, is driven by the magnetic pressure gradient. This 
has been proposed elsewhere^^' — , but a quantitative comparison has so far been lacking. 
We have shown with a PIC simulation that the force imposed on an electron by the time- 
averaged electrostatic field Ex{x) matches the Eb{x), which results from the time-averaged 
magnetic pressure gradient force. The Ex{x,t) is, however, not time-stationary, which can 
be explained as follows. 

The FI accelerates through the magnetic pressure gradient force the electrons and a 
current Jx{x, t) builds up. This current results in the ID geometry with V x S = through 
the normalized equation Jx{x,t) = —dtEx{x,t) in a growing Ex{x,t). The initial conditions 
are Ex{x,t = 0) = and Jx{x,t = 0) = 0. Any oscillatory solution for Jx{x,t) and Ex{x,t) 
implies through Jx{x, t) = —dtEx{x, t) that Jx and Ex cannot simultaneously oscillate in time 
around their initial values. The Ex{x,t) oscillates instead around its time-average, which is 
the background field Eb{x). The oscillation amplitude of Ex is approximately Eb{x). The 
superposed oscillatory and background field thus oscillates between Ex{x,tc) = at certain 
times tc, fulfilling the initial condition at = 0, and a maximum Ex{x,tc) = 2Eb{x). 

We have confirmed previous suggestions, that the electric field force is comparable to 
the magnetic field forced' —. We have used the time-averaged electric and magnetic forces 
to estimate their effects quantitatively. The electric field repels electrons at the filament 
centres and attracts them if they are farther away, which permits filaments to overlap and 
limits their peak density. 
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